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ABSTRACT 


The purpose of this research was to demonstrate the feasibility of porting a 
large mainframe-dependent scientific FORTRAN program, specifically the 
Numerical Electromagnetics Code (NEC-3) to a 32-bit personal computer. Two 
systems, an AST Premium 386/33 with both Intel 80387 and Weitek w3167 math co- 
processors and a Definicon DSI-780 using a Motorola 68020 CPU and 68881 math 
co-processor, were used with several 32-bit FORTRAN 77 compilers. 

Results show that when NEC-3 was promoted to full double precision, 
complete accuracy was maintained while suffering only a 12% increase in execution 
time over single precision. Testing also revealed that the double precision Weitek 
version is 30% faster than the 80387 version. Some small inaccuracies remain; 
however, the same results were obtained by both the Naval Postgraduate School 
mainframe’s new IBM VS2 FORTRAN 77 compiler and the personal computer 
FORTRAN 77 compilers. This indicates that the "bug" is in the NEC-3 code, vice 


being a hardware or FORTRAN 77 compiler problem. 
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I, INTRODUCTION 


A. PURPOSE 

The Numerical Electromagnetics Code version-3 (NEC-3) advanced antenna 
modeling program was first implemented on the Naval Postgraduate School IBM 
3033AP mainframe computer in 1983. Since then, NEC-3 has proven to ‘e an 
invaluable tool in antenna design, particularly in the high interest area of buried 
antennas. The use of a personal computer version of NEC-3 will greatly benefit 
researchers by freeing them from the constraints of the mainframe computer. 

The principle objective of the research was the evaluation of the performance 
of a large mainframe-dependent FORTRAN program, specifically NEC-3, which had 
been ported to a fast 32-bit personal computer. Recent advances in microcomputer 
architecture, semiconductor memory, and hard disk drive performance now make it 
practical to implement such large scale scientific FORTRAN programs on a pe’sonal 
computer. Technically, the mainframe is still faster than most 32-bit personal 
computers, although this speed advantage is lost if more than a few users are on the 
system. The relatively inexpensive cost, reliability, portability, and user friendliness 
have made the 32-bit personal computer the ideal computer system for scientific 


programming and computation. 





B. BACKGROUND 


1. Mainframe Computer Versions 

The NEC-3 is a user-oriented FORTRAN program for the analysis of the 
electromagnetic response of antennas and other metal structures. NEC-3 is the latest 
version of a series of electromagnetic FORTRAN programs, each of which improved 
upon its predecessors. The first program was BRACT, developed by MBAssociates 
in San Ramon, California in 1970, under the funding of the Air Force Space and 
Missiles Systems Organization. BRACT was specialized to scattering by arbitrary 
thin-wire configurations. In 1974, the Antenna Modeling Program (AMP) was 
developed by MBAssociates with funding from the Naval Research Laboratory, 
Naval Ships Engineering Center, U.S. Army ECOM/Communications Systems, U.S. 
Army Strategic Communications Command, and the Rome Air Development Center. 
AMP improved upon BRACT in several aspects by adding the capability of modeling 
a structure over a ground plane and allowing the use of file storage, thereby 
increasing the maximum structure size that could be modeled. AMP-2, introduced 
in 1975, included a simplified approximation for large interaction distances which 
helped reduce run time for large structures. [Ref. 1:p. 1] 

NEC-1 was developed at the Lawrence Livermore Nationa] Laboratory 
(LLNL) under the sponsorship of the Nava] Ocean Systems Center (NOSC) in 1977. 
NEC-1 added to AMP2 a more accurate current expansion along wires and at wire 


junctions, and an option in the wire modeling technique for greater accuracy for 








thick wires. A new model for a voltage source was included along with several! other 
modifications made for increased accuracy and efficiency. NEC-2 was introduced 
in 1980, with the added capability of performing the Numerical Green’s Function for 
partitioned-matrix solutions and a treatment for lossy grounds that is accurate for 
antennas very close to the ground surface. In 1983, the current version, NEC-3, was 
released. This version has the added capability of modeling buried antennas and 


wires penetrating a lossy media. [Ref. 1:p. 2] 


2. Personal Computer Versions 

The personal computer’s relatively low cost and increasing capabilities 
have made it the focus of several efforts to implement a working version of NEC 
for a personal computer. 

The NOSC developed a Mini-Numerical Electromagnetics Code 
(MININEC) in 1982. This program was written in BASIC and included many, but 
not all of the capabilities of the full NEC program. MININEC was also limited to 
modeling only electrically small wire structures by the amount of memory available 
(640 KB) in the personal computer. [Ref. 2:p. 1] 

Stephan Lamont of the Naval Postgraduate School (NPS) was the first 
known researcher to implement the full version of the NEC-3 program on small 
computer systems in 1986. The three systems that he used were: 


e« An IBM PC-AT, as an example of a commonly available "off the shelf" 
microcomputer. 


e An IBM PC-AT/370, as a middle-to-high end micro-mainframe 

e An IBM RT PC, as an example of a high-performance, high-function 
workstation embodying a novel architecture and with the capabilities outside 
the range of personal computers. [Ref. 3:pp. 1-2] 

The results of his research were somewhat promising and showed that for 
the first time it was possible to run the full NEC-3 program on a small "personal" 
system. Both the IBM RT PC and IBM PC AT/370 executed without any problems; 
the IBM PC AT, however, suffered from several problems which are summarized 
as follows: 

¢ The IBM Professional FORTRAN did not support the DOUBLE PRECISION 
COMPLEX data type, thus requiring double precision complex algorithms to 
be coded in two subroutines to provide the required accuracy. 

e When processing an out-of-core solution, a record length error frequently 
occurred, requiring the addition of a run-time parameter specifying the 
maximum record length allowable. 

e Out-of-core solutions did not produce correct results. 

¢ The speed of the IBM PC-AT was between 2 and 10 times slower than that of 
the IBM RT PC or IBM PC-AT/370. [Ref. 3:pp. 3-6] 

Captain Timothy O’Hara, U.S. Army, implemented the full NEC-3 
program on three 32-bit PC systems in 1988: 

¢ a COMPAQ 386/20 with both Intel 80387 and Weitek w1167 math co- 
processors, 


e an IBM RT PC with a National Semiconductor (NS) 32091 floating point 
accelerator, and 








e a Definicon DSI-780 co-processor board with a Motorola 68020/20 and a 
Motorola 68881 math co-processor. 

All three systems had large in-core memories, with 4 Megabytes each. The 32-bit 
MicroWay FORTRAN compiler used on the COMPAQ 386/20 implemented 
extended memory via a DOS Extender made by Phar Lap Tools. [Ref. 4:pp. 3-4] 
The results of O’Hara’s research showed that although both the Definicon 

and IBM RT PC systems generated accurate data for simple problems, their 
execution times were not fast enough to justify their use for large-scale problems. 
The COMPAQ 386/20, however, was fast enough but had several very serious 
accuracy problems. When the Intel 387 math chip was used, the antenna input 
impedances became stochastic because of a COMMON block alignment problem. 
When the Weitek w1167 was used, the results became more inaccurate as the 
number of unknowns in the problem increased. The Weitek inaccuracy problem was 
attributed to the small single precision mantissa of the Weitek chip, 23 bits versus 53 


bits in double precision. [Ref. 4:pp. 33-34] 


C. CONTENTS 
Chapter II contains a detailed description of the AST Premium 386/33 personal 
computer, the Definicon DSI-780 co-processor board, and the FORTRAN compilers 


used by each system. 











Chapter II] contains a description of past and current FORTRAN language 
standards, their impact upon the NEC-3 program, and the specific changes that were 
made to NEC-3 during the research. 

Chapter IV contains a description of the testing of the latest personal computer 
version of NEC-3 and the results. 


Chapter V contains conclusions and recommendations. 








II. PROGRAMMING ENVIRONMENT 


Two separate 32-bit computer systems were utilized for this research, an AST 
Research Premium 386/33 MS-DOS compatible system and a Definicon DSI-780 co- 
processor board. Three different FORTRAN compilers were used on the MS-DOS 
system in order to provide a ready means of completing the project in the event that 
any one compiler was unsuccessful in running the NEC-3 program. Only one 
FORTRAN compiler was available for the Definicon DSI-780 co-processor board, 


and the most recent version of that compiler was used. 


A. HARDWARE 


1. AST Premium 386/33 
The AST Premium 386/33 System is a very high performance MS-DOS 
system using an Intel 80386 microprocessor running at 33 MHz. The system has 4 
MB of 32-bit memory on the motherboard, 8 MB of 32-bit memory installed on an 
AST memory expansion board, and 32 KB of cache memory. Both Intel 80387 
running at 33 MHz and Weitek w3167 floating point co-processors are installed on 
the motherboard. Hard disk storage is provided by a fast EDSI 100 MB hard disk, 


and the Microsoft SMARTDRV.SYS disk cache routine is used to improve disk 


performance throughput. The Norton Utilities computing index (C1) for this system 


relative to the IBM/XT is 38.8. (Ref. 5:p. xi] 


2. Definicon DSI-780 Co-processor Board 
The Definicon DSI-780 is a 32-bit microprocessor, running at 20 MHz, a 
32-bit data bus, an MC 68881 floating point math co-processor, and 4 MB of 32 bit 
dynamic memory (RAM). The board interfaces the MS-DOS operating system of 


PC-AT computers via a 64 KB window located at DOOOH. [Ref. 6:pp. 1-3} 


B. FORTRAN COMPILERS 

The three MS-DOS FORTRAN compilers selected for this research were 
chosen because they support the use of extended memory and the DOUBLE 
PRECISION COMPLEX data type. The extended memory feature was a 
requirement because of the large program size and in-core data storage requirements 
of the NEC-3 program. The DOUBLE PRECISION COMPLEX feature was 
required to produce a full double precision version of NEC-3. This was desirable 
because the IBM mainframe version was promoted to full double precision at 
compile time, and the single precision version run by O’Hara had some precision 
related accuracy problems. 

The only FORTRAN compiler available for the Definicon DSI-780 co- 
processor board was the SVS FORTRAN compiler. The newest version, v.2.8, was 


used as it now supports the DOUBLE PRECISION COMPLEX data type. Memory 








access on the Definicon DSI-780 was not a problem since it was not affected by the 


MS-DOS memory limit. 


1. Lahey F77L-EM/32 FORTRAN Compiler (v.2.1) 

The F77L-EM/32 32-bit FORTRAN compiler v.2.1 is produced by Lahey 
Computer Systems, Inc. of Incline Village, Nevada. This compiler fully supports the 
Intel 80386 microprocessor, Intel 80387 math co-processor, and the Weitek w3167 
math co-processor. Extended memory support is provided via a transparent DOS 
Extender which requires no special action by the program to use extended memory. 
[Ref. 7:p. i] 

The F77L compiler is a true native code compiler that compiles directly 
to 80386 machine code without going through an assembly language first. The 
compiler code is linked into an executable file with an (.EXP) extension by the 
Lahey linker (L32). The executable file requires the Lahey loader (UP) in order to 
execute properly in extended memory. Appendix A contains a complete listing of 
all switches and options available with the Lahey F77L FORTRAN compiler. Of 
special note is the "/Q1" switch. This switch forces the compiler to limit the math 
co-processor stack to eight registers, which is the physical limit of the math co- 
processor. The default case is for this switch to be turned off, which allows the math 
co-processor stack to grow beyond eight registers via an exception trap to the 
runtime library scheme. When this happens while executing NEC-3, the results 


were completely inaccurate and unusable. It is unclear why an unlimited math 




















co-processor stack would be useful to anyone. The documentation only indicated that 
the switch rarely causes any problems, which in this case includes the AST 386/33. 


NEC-3 executed properly when the switch was turned on. [Ref. 7:pp. A1-A8] 


2. MicroWay NDP FORTRAN-386 Compiler (v.1.4E) 

The 32-bit NDP FORTRAN-386 compiler v.14E is produced by 
MicroWay, Inc.,of Kingston, Massachusetts. Full support is provided for the Intel 
80386 microprocessor, Intel 80287 and 80387 math co-processors, and the Weitek 
w3167 math co-processor. Full access to extended memory is provided with the Phar 
Lap Tools DOS Extender. The DOS extender, however, is not part of the NDP 
FORTRAN-386 package and must be purchased separately in order to produce 
executable files. The DOS extender is transparent to the programmer, and extended 
memory is fully available to the program at execution time. [Ref. 8:pp. 1.1-1.2] 

NDP FORTRAN-386 compiles to an intermediate assembly language code, 
which is then assembled by the Phar Lap Tools assembler, ASM386, into object 
code. The Phar Lap Tools linker, LINK386, produces executable (.EXP) file from 
the object file. The Phar Lap Tools loader, LOAD386, is required to load the 
executable file into extended memory. Appendix B contains a complete list of all 


compiler switches and options. [Ref. 8:pp. 1.5-1.12] 


10 








3. SVS FORTRAN-386 Compiler (v.2.8) 

The SVS FORTRAN-386 compiler v.2.8 was distributed by Science 
Applications International Corporation (SAIC) of Los Altos, California. The 
compiler is produced by Silicon Valley Software (SVS) of Cupertino, California. 
The linker is produced by Phar Lap Tools and the loading executive manager is 
provided by Intelligent Graphics Corporation (IGC) of Santa Clara, California. Full 
support for extended memory is provided by the VM/RUN loading execution 
manager. All DOS extended memory is available to the program at execution time. 
[Ref. 9:p. 1.1] 

SVS FORTRAN-386 is a three pass compiler which uses three separate 
programs to compile the source file into an object file (.OBJ). The object file is 
then linked into an executable file by the Phar Lap Tools Linker (LINK386). The 
loading execution manager, VM/RUN, is required to execute the file in extended 
memory. Appendix C contains a list of all compiler switches and options. [Ref. 


9:pp. 3.1-3.2] 


4. SVS Definicon FORTRAN Compiler (v.2.8) 

The SVS Definicon FORTRAN compiler (v.2.8) is produced by SVS of 
Cupertino, California. The linker is also made by SVS, and the loader is provided 
by the co-processor board manufacturer, Definicon Systems, Inc. The 32-bit 
FORTRAN compiler fully supports the Motorola MC 68020 Microprocessor, and the 


Motorola MC 68881 floating point co-processor. [Ref. 10:pp. E.1-E.3] 








The SVS Definicon FORTRAN compiler is a two-step compiler where the 
source file (.FOR) is first compiled to an intermediate step before being compiled 
to an object file (OBJ). This is in contrast to the SVS MS-DOS FORTRAN 
compiler which has two intermediate steps. The resulting object file is linked into 
an executable file by the linker. A special nuance of the Definicon system is that 
the loader must be used to move al! of the 68020 executable file (compiler, linker, 
program) from the MS-DOS host computer to the Definicon DSI-780 board. The 
speed of the loader is directly related to the speed of the host system. Appendix C 


contains a list of all compiler switches and options. [Ref. 10:pp. E.1-E.3] 
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Ill. PROGRAM SOURCE CODE 


A. THE FORTRAN PROGRAMMING LANGUAGE 

The first American National Standard for the FORTRAN programming 
language was established in 1966 by the American National Standards Institute 
(ANSI) and was known as FORTRAN 66. However, most FORTRAN compilers 
that were based on this standard went by the name FORTRAN IV. The standard 
permitted the use of "compatible extensions,” and by the mid 1970’s it became 
apparent that many of the more commonly used extensions to FORTRAN 66 could 
be incorporated into a new FORTRAN language standard. In 1978, ANSI formally 
upgraded the FORTRAN language standard to what is commonly referred to as 
FORTRAN 77. Although the standard is voluntary, all modern FORTRAN 
compilers fully support the FORTRAN 77 standard. [Ref. 11:p. 427] 

Most older FORTRAN programs, in particular NEC-3, were written before 
1978. Because of this, it is important to understand the differences between 
FORTRAN 66 and FORTRAN 77, and how these differences will affect program 


portability. 








1. Major Differences Between FORTRAN 77 and FORTRAN 66 


a. Structured Branching Statements 
The IF-THEN-ELSE-END IF control structure has been added to the 
language. Unlimited nesting of these control structures is allowed. Nesting within 


DO loops is also permitted. [Ref. 11:p. 480] 


b. Character Data Type 
A new data type, consisting of character strings of fixed declared 
length, has been added to the language, including character constants, character 
variables, and arrays of character data. Concatenation and designation of character 
substring is also supported. The Hollerith data type of ANSI x3.9-1966 has been 


deleted. [Ref. 11:p. 481] 


c. DO Loop Changes 
A DO statement specifying a terminal parameter whose value is less 
than that of the initial parameter is no longer prohibited, under such conditions the 
loop will be skipped and not executed. FORTRAN 66 specified that all loops must 
execute at least once. Negative increments are also permitted. Transfer of control 


into a DO loop is now prohibited. [Ref. 11:p. 481] 
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d. List Directed Input and Output 
A form of input and output is provided which does not require an 


explicit format specification. [Ref. 11:p. 481] 


e. Expressions 
An arithmetic expression may include subexpressions of more than 
one type. The operands are first converted to the data type of the result. DO loop 
parameters can be any expression of integer, real, or double precision type. [Ref. 


11i:p. 481] 


f. Compile-time Constants 
The PARAMETER statement has been added, which declares the 
value corresponding to the symbolic name of a constant. Such a name may be used 
in an expression, DATA statement, or in following PARAMETER statements. [Ref. 


11:p. 481) 


g- IMPLICIT Type Declaration 
An IMPLICIT statement may be used to modify the default typing 


by the compiler for variable and array names beginning with certain letters. [Ref. 


11:p. 482] 














h. Generic Intrinsic Functions 
Most intrinsic functions such as SIN and COS now have a generic 
name, so that special names for complex and double precision (i.e., CSIN, DSIN) 


need not be used. [Ref. 11:p. 482] 


i. Subprogram Reference 
Subroutines and functions may contain ENTRY statements, and the 


alternate return may be used in subroutines. [Ref. 11:p. 482] 


j. Array Bounds 
Array declarations may include both upper and lower dimensions. 
Arrays may have up to seven dimensions. The last dimension of an array passed to 
a subroutine may be dimensioned inside the subroutine with an asterisk (*). 
FORTRAN 66 specified that the last array dimension of an assumed array be 


dimensioned with a one (1). (i.e., A(10,1) vs. A(10,*)) (Ref. 11:p. 482] 


k. Computed GOTO Default 
If the control expression of a computed GOTO is out of range, 
execution continues with the statement following the computed GOTO. [Ref. 11:p. 


482] 
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1. INPUT/OUTPUT Statements 


Several new features were added to the FORTRAN 77 standard to 


improve INPUT/OUTPUT. [Ref. 11:p. 482] 


An OUTPUT list may now contain constants and expressions. 


INPUT/OUTPUT statements may now contain a character string to be used 
as a format specification. 


END OF FILE and ERROR condition control. 
Tab format edit descriptors have been added. 
Direct access file support has been added. 
Character arrays maybe used as an internal file. 


OPEN, CLOSE and INQUIRE statements added. 


m. SAVE Statement 


The SAVE statement is now part of the standard. This statement 


instructs the compiler to retain the value of specified local variables between calls 


to the subprogram. [Ref. 11:p. 482] 


n. FORTRAN Character Set/Comment Lines 


The apostrophe (’) and colon (:) are now part of the standard 


FORTRAN character set. Both the asterisk (*) and letter C in column 1 designate 


a comment line. [Ref. 11:p. 482] 











2. Portability of FORTRAN Programs 

FORTRAN is one of the few universal programming languages. For 
several decades it has been the primary computer language for scientific, numerical, 
and technical applications. Only now are some inroads into its position being made 
as simpler, better structured languages such as PASCAL, C, and ADA become more 
popular for numerical applications. FORTRAN’s universality is based on the wide 
availability of compilers for virtually all sizes of computers from the super-computer 
to the microcomputer. While this has helped the exchange of FORTRAN programs, 
there have been portability problems. [{Ref. 12:p. 134] 

Many compiler manufacturers have introduced extensions to the language 
standard. These extensions fall into two broad categories, those which are true 
extensions of the power of the language such as new types of control structures and 
those which allow access to special hardware features of the system the compiler is 
supporting. Both types of extensions work against program portability, the ability to 
transfer programs from one computer to another with no source code modifications 
and to obtain the same results and level of precision. [Ref. 12:p. 134] 

A second barrier to program portability is the FORTRAN standard itself, 
since not all parts of the language specifications are completely defined. For 
example, the depth to which DO loops can be nested is compiler-dependent. Such 
differences between compilers can lead to difficulties when moving programs to new 


computer systems. [Ref. 12:p. 135] 
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The best method of guaranteeing program portability is to strictly adhere 
to the standard when the program is created. This rarely occurs, however, since 
usually only one compiler is available when the program is written and the 
peculiarities of the compiler always seem to find their way into the program. There 
are six categories of the FORTRAN standard where most portability problems occur, 
language elements, expressions and assignments, contro] statements, specification 


Statements, program units, and INPUT/OUTPUT. (Ref. 12:p. 135] 


a. Language Elements 


(1) Character Set. The standard defines upper case letters only, 
lower case should not be used except in comment lines. In addition, there are 
thirteen special characters: 

blank,.’: = +-*/()$ 
and no other non-alphanumeric characters should be used. [Ref. 12:p. 136] 
(2) Source Form. The standard defines a rigid source form, columns 


1:5 for the label, column 6 for the continuation mark, columns 7:72 for the 


statement. This format should be respected. [Ref. 12:p. 136] 


(3) Continuation Lines. The standard prohibits the presence of any 


non-blank characters in columns 1:5 of a continuation line; not all compilers enforce 


this, but the field should be left blank anyway. (Ref. 12:p. 136] 











(4) Multiple Statement. Some compilers allow more than one 
statement to be written on a line, usually separated by the non-standard semi-colon 


(;). Not only is this non-portable, it is bad programming style. [Ref. 12:p. 137] 


(5) In-line Comments. Some compilers allow comments to be added 
to statements on the same line following the non-standard exclamation mark (!). 


This is a UNIX style extension which should be avoided. [Ref. 12:p. 137] 


(6) Long Names. Some compilers allow names to have more than 
six characters. While this is an improvement in style, it is also highly non-portable 


and very difficult to modify later. [Ref. 12:p. 137] 


(7) Currency Symbol. The standard defines no particular role for 
the currency symbol ($). Some compilers, however allow it to be used as a character 


in names. This again is non-portable and should be avoided. (Ref. 12:p. 137] 


(8) Logical Constants. The logical constants are .TRUE. and 
-FALSE.. Some compilers allow abbreviations such as .T. and .F. which should be 


avoided. [Ref. 12:p. 137] 


(9) Character Strings. Character strings are delimited by 
apostrophes (’). No other character, in particular the non-standard double quotes ("), 


should ever be used. [Ref. 12:p. 138] 
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b. Expressions and Assignments 
This section deals with some potential problems with expressions and 


assignments, mostly related to different hardware implementations of FORTRAN. 


(1) Mixed Mode Relational Expressions. Although mixed mode 
relational expressions such as 
IF (A/B .NE. I-J) THEN 
are permitted by the standard, they should be avoided since real expressions will 


evaluate slightly differently on different machines. [Ref. 12:p. 138] 


(2) Tests of Accuracy. It is also important to avoid tests of equality 
involving floating point operands, as hardware dependent rounding errors may cause 


them to be evaluated incorrectly. (Ref. 12:p. 138] 


(3) Subscript Expressions. Some FORTRAN compilers allow the 
use of non-integer subscript expressions. The value of the expression is rounded 
according to the standard rules. Because rounding errors can cause different results 
on different machines, inherently real expressions should be rounded explicitly to the 
nearest integer value using the NINT intrinsic function, as in 

A (NINT (P/Q)) = 0. 


[Ref. 12:p. 139] 
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(4) Number of Subscripts. The number of subscripts in an array 
reference should be identical to the number of dimensions of the array as in 
INTEGER | J(10,10) 
j (3,4) = 0 
and not: 
J (33) = 0. 


{Ref. 12:p. 139] 


c. Control Statements 
This section covers several areas where control statements can 


contribute to the non-portability of a program, and which should always be avoided. 


(1) DO-loop Parameters. FORTRAN 77 permits the use of real and 
double precision quantities as DO-loop indices and parameters, both in evaluating 
the number of interactions and the successive values of the loop index. [Ref. 12:p. 


140] 


(2) Arrays in Nested Loops. Nested loops, when used, should be 
arranged so that the subscripted array’s first element depends on the innermost loop 


and its last on the outermost loop, such as: 


DO1 J = 1, 10 
DO2 I = 1, 10 
A (I,J) = B (IJ) + C (I,J) 
2 CONTINUE 
1 CONTINUE 
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When loops are arranged in this manner, unexpected and dramatic losses in 
efficiency connected with the way words are stored in main memory can be avoided. 


[Ref. 12:p. 140] 


(3) Loop Constructs. Many compilers now allow the use of 
PASCAL-type loop constructs, such as DO... WHILE and DO. . UNTIL. These are 


non-portable and should never be used. [Ref. 12:p. 140] 


(4) Branches into Control Constructs. Some compilers allow control 
to be passed into an IF-THEN block. This is bad style as well as being non-portable. 


[Ref. 12:p. 141] 


d. Specification Statements 


The section deals with the rules concerning specification statements. 


(1) Data Types. The standard defined six data types, INTEGER, 
REAL, DOUBLE PRECISION, COMPLEX, LOGICAL, and CHARACTER. No 
other data types should be used, nor should non-standard declarations, such as 
REAL *4, or REAL*8 be used when a standard declaration has the same meaning. 


{Ref. 12:p. 141] 


23 














(2) EQUIVALENCE Statements. The EQUIVALENCE statement 
should never be used for variables of different types, since internal number 
representations differ from computer to computer, resulting in non-portable code. 


{Ref. 12:p. 141] 


(3) PARAMETER Statements. Some compilers allow intrinsic 
function calls to be placed in constant expressions inside PARAMETER statements. 


This is not part of the standard and should be avoided. (Ref. 12:p. 141] 


(4) DATA Statements. FORTRAN programs should never rely on 
the default data initialization of a specific machine/operating system combination. 
(for instance, the setting of all variables not explicitly defined by a DATA statement 
to zero). All variables should be explicitly set in either an assignment statement or 


DATA statement. [Ref. 12:p. 141] 


e. Program Units 
This section cover program units and the specification statements 


which are relevant to them. 


(1) Main Program Header Line. All complete FORTRAN programs 
must contain a main program. The main program header line is optional, but if used 
it should not contain a parameter list. The parameter list was used in FORTRAN 
66 as a means of defining the I/O files. All such files should be explicitly defined 


by the OPEN statement. [Ref. 12:p. 142] 
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(2) Intrinsic Functions. Most compilers offer additional intrinsic 
functions beyond those defined by the standard. These additional functions should 
never be used, since even if the same function were to be available on another 
system, the syntax would most likely be different. The most commonly available 


additional function is the TIME/DATE function. [Ref. 12:p. 142] 


(3) Character-to-Integer Conversion. The values returned by the 
intrinsic functions CHAR and ICHAR are processor-dependent, and no portable 


program should depend on any specific values. [Ref. 12:p. 142] 


(4) External Functions. Although not required to do so by the 
standard, all external functions should be explicitly declared in an EXTERNAL 
statement to avoid surprises. This is because the standard allows the compiler to 


assign any known function to an undeclared reference. [Ref. 12:p. 142] 


(S) User-Defined Functions. User-defined functions should not alter 
the value of their arguments, or COMMON variables, nor perform 1/O. The order 
in which any of these undesirable actions is performed is compiler dependent and 


leads to non-portable code. [Ref. 12:p. 123] 


(6) COMMON Variables as Arguments. COMMON block variables 
should not be passed as arguments when the COMMON block is referenced in both 
the calling and called subroutine. It is especially hazardous to pass COMMON block 


variables when they are modified inside the called subprogram. [Ref. 12:p. 143] 
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(7) CHARACTER Variables in COMMON Blocks. The standard 
does not define any mapping of CHARACTER variables onto computer words as 
it does for other data types. Although permitted by many compilers, the mixing of 
CHARACTER variables with other data types in common blocks is highly non- 
portable since the size of a character storage unit varies from computer to computer. 


[Ref. 12:p. 143} 


(8) Double-word Data Types in COMMON Blocks. Some compilers 
require that large data types, double precision and complex, should be positioned in 
memory so as to begin on an odd word. Since this is a grey area in the standard, all 
such variables should be positioned at the beginning of the common block so they 


will always be guaranteed of beginning on an odd word. [Ref. 12:p. 143] 


(9) Names. Although allowed by some compilers, no COMMON 


block name should be the same as the name of any subprogram. [Ref. 12:p. 144] 


(10) Order of Statements. While many compilers are lax in enforcing 
the order of statements as defined by the standard, all portable programs should 


strictly follow the correct order of statements. [Ref. 12:p. 144] 


(11) SAVE Statement. The SAVE statement should be used as 
appropriate, as it might be required by some compilers on which the program will 


Jater run. (Ref. 12:p. 144] 
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(12) BLOCK DATA. When DATA statements are used for the 
initialization of COMMON variables, a BLOCK DATA subroutine should be used. 


[Ref. 12:p. 144] 


f. Input/Output 
This section deals with several points about I/O which can affect 


program portability. 


(1) PRINT and READ statements. PRINT and READ statements 
without a unit specifier cause I/O to be directed to or from compiler-defined units, 
and, as such, they should not be used in portable code. WRITE and READ 
statements should use variables or symbolic constants for unit specifiers, since they 


can be easily modified if necessary. [Ref. 12:p. 144] 


(2) Non-standard 1/O. Many compilers have non-standard 1/O 
statements which should never be used in portable programs. Non-standard 


parameters in standard I/O statements should also be avoided. (Ref. 12:p. 145] 


(3) Parameter Lists. Both READ and WRITE statements have 
either positional or keyword parameters. Both styles should never be mixed in one 


statement. [Ref. 12:p. 145] 
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(4) I/O Unit Numbers. I/O unit numbers should never be explicit 
constants, or an asterisk {*), but rather symbolic constants or variables which can be 


easily changed when moving code from one system to another. [Ref. 12:p. 145] 


(5S) Format Specifications. Format Specifications may be stored in 
arrays. To avoid portability problems the arrays should always be of type 
CHARACTER. Additionally, format specifiers should always be separated by a 


comma (,) to ensure portability. [Ref. 12:p. 145] 


(6) IOSTAT Parameter. The values returned by the IOSTAT 
parameters in I/O statements are compiler dependent and should never be explicitly 


enumerated; symbolic constants should be used instead. (Ref. 12:p. 145] 


B. OVERVIEW OF THE NEC-3 FORTRAN PROGRAM 


1. Theory 
The NEC-3 is a highly capable, and versatile tool for the analysis of the 
electromagnetic response of antennas and other metal structures. The program 
computes directly the numerical solution of integral equations for currents induced 
on the structure by sources or incident fields, thus avoiding many of the simplifying 
assumptions required by other solution methods. [Ref. 1:p. 1] 
The electric-field integral equation (EFIE) and the magnetic-field integral 


equation (MFIE) are both used to model the electromagnetic response of general 
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structures. Each equation has different advantages for a particular structure type. 
The EFIE is best suited for thin-wire structures of small or vanishing conductor 
volume. The MFIE, which does not work for thin-wire structures is most useful 
when modeling voluminous closed surface structures, particularly those with large 
smooth surfaces. The EFIE and MFIE are coupled for structures which contain both 
surfaces and wires. [Ref. 1:p. 3] 

The integral equations are solved numerically using a form of the Method 
of Moments. The resulting current matrix equation is then solved by Gauss 
elimination and LU decomposition. The computation of the matrix elements and the 
solution of the matrix equation are by far the two most time-consuming steps in 
computing the response of a structure, usually accounting for 90% or more of the 
total computation time. [Ref. 1:p. 30-32] 

The order of the matrix equation (N) is directly affected by the number 
of wire segments (N,) and surface patches (N,) in the structure being modeled. [Ref. 
i:p. 30] 

N = N, + 2N, (IIl-1) 

When the structure is modeled near ground, penetrating the ground, or buried, 

the SOMNTX FORTRAN program is used to generate a data file containing the 
relevant Sommerfeld integral interpolation tables. By using a separate program to 
generate these interpolation tables,the computation time for ground related problems 


is greatly reduced. [Ref. 1:p. 37] 
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2. Program Organization 

The NEC-3 program consists of over 10,000 lines of FORTRAN code, 89 
subroutines and user-defined functions, and over 1000 variables. The program is 
entirely self-contained, i.e., no libraries are required. The SOMNTX program is 
somewhat more modest, being only 2000 lines long with only 25 subroutines and 
user-defined functions. Both programs were written largely in FORTRAN 66 and 
use archaic programming elements (GOTO, computed goto, assign, equivalence) 
extensively. Very little use is made of the advanced FORTRAN 77 control 
structures such as the IF-THEN-ELSE block. As a result, program flow is very 


difficult to follow. 


C. MODIFICATIONS INCORPORATED INTO THE NEC-3 FORTRAN 
PROGRAM 


The following are the modifications that have been incorporated into the NEC- 
3 and SOMNTX FORTRAN source code as part of this research. Most of the 


changes are related to the topics covered in the program portability section. 


1. NEC-3 Modifications 


a. SAVE Statements 
Gerry J. Burke of the Lawrence Livermore Laboratory discovered 
that SAVE statements were required for all subprograms when he ported the VAX 


version to the MACINTOSH II. The justification for this is that in FORTRAN 66, 
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most compilers save all local variables automatically. By adding the SAVE 


statement to the NEC-3 the stochastic problem reported by Tim O’Hara disappeared. 


b. Common Block Alignment 
When the UNIX version of the MicroWay NDP FORTRAN compiler 
was used, it was discovered that a common block misalignment problem existed. 
The complex variable T1 was not defined as type COMPLEX in those subroutines 
where the variable was not referenced. However, since the common block GND 
was a part of those subroutines, it is mandatory that the variable T1 be explicitly 
declared as type COMPLEX in all subroutines where the common block GND is 


used. 


c. Full Double Precision 
The NEC-3 program was promoted to full double precision based on 
the recommendations of Tim O’Hara’s research. This was done by placing the 
IMPLICIT REAL *8 (A-H, O-Z) statement at the beginning of all subprograms, and 


making the following intrinsic function replacements: 


REAL —- DBLE 
FLOAT ~ DBLE 
COMPLEX -~+ COMPLEX*16 


CMPLX -~ DCMPLX 
AIMAG - DIMAG 


CONJG DCONJG 


4 
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d. Generic Functions 
As mentioned earlier, one of the improvements of FORTRAN 77 
over FORTRAN 66 was the introduction of generic functions. To make the 
conversion of NEC-3 to double precision easier, all intrinsic functions except DBLE, 
DCMPLX, DIMAG, AND DCONJG were converted to their generic functions. For 


example, CABS became ABS. 


e. FORTRAN 77 Standard Array Passing 
FORTRAN 77 specifies that all assumed size dummy arrays use an 
asterisk, e.g., A(*), rather than the FORTRAN 66 method of dimensioning assumed 
size arrays with the number one (1), e.g., A(1). Since the old method was used 
extensively in NEC-3 and is not documented by any of the compilers used for this 
research, the author felt it prudent to convert all assumed size arrays to the 


FORTRAN 77 standard. 


f. Output Statements 
As covered previously, output statements without a unit specifier 
direct output to the default device, which is compiler-specific. In NEC-3, the PRINT 
statement sends the output to unit 6 when using the MicroWay compiler and to the 
screen when using the Lahey compiler. To ensure portability all PRINT statements 


were changed to WRITE statements. 
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g- External Functions 
All external functions used in a program unit are now first declared 


in an EXTERNAL statement. 


h. Subroutine Names 
The names of two subroutines, INPUT and READ, conflicted with 
the names of routines in the SVS compiler’s library. The solution was to rename the 


subroutines as INPUT1 and READ1, and also modify all the CALL statements. 


i. INTEGER*4, REAL*8, and COMPLEX *16 

All compilers make undeclared integer variables to be of type 
INTEGER*4, however, because all MS-DOS compilers now support the 
INTEGER*2 data type, the INTEGER statement is incorrectly interpreted by the 
compiler to be of type INTEGER*2 vice INTEGER*4. Severe problems occur 
whenever 4-byte integers are passed to 2-byte integers and vice versa. REAL*8 and 
COMPLEX* 16 were used in lieu of DOUBLE PRECISION REAL and DOUBLE 
PRECISION COMPLEX, respectively, since all modern FORTRAN 77 compilers 


support these extensions. 


j. Arrays of Mixed Data Types 
In an effort to save space, the NEC-3 program reuses arrays 
whenever possible. Most of the time, real variables replace real variables, and no 


harm is done. In the DATA and DATAJ common blocks several integer arrays are 
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used to hold both integer and real variables. When single precision real variables 
are used, both real and integer variables are 4 bytes long, and can be used in the 
same array via an EQUIVALENCE statement without causing problems. However, 
this scheme does not work when double precision real variables are used. 

In order to resolve conflicts, the REAL*8 arrays, T2X(1000), 
T2Y(1000), T2Z(1000), have been added to the COMMON block DATA. The 


following DIMENSION and EQUIVALENCE statements have been deleted: 


e DIMENSION T2X(1000), T2¥(1000), T2Z( 1000) 
e EQUIVALENCE (T2X(1), ICON1(1)), (T2Y(1), ICON2(1)), (72Z(1), 
ITAG(1)) 
The affected subroutines are: MAIN, ARC, CABC, CMNGF, CMSET, CMSS, 
CMSW, CMWS, CMWW, CONECT, DATAGN, ETMNS, FFLD, FFLDS, GFIL, 
GFLD, GFOUT, ISEGNO, LOAD, MOVE, NEFLD, NETWK, NFPAT, NHFLD, 
PATCH, QDSRC, RDPAT, REFLC, SBF, TBF, TRIO, and WIRE. 


k. EQUIVALENCE Statements with Mixed Data Types 
Like arrays, the mixing of data types of different lengths can cause 
severe problems. Two INTEGER variables in the DATAJ COMMON block were 
equivalenced to REAL*8 variables part of the time. The following changes were 


made: 








e Add the REAL*8 variables T2YJ and T2ZJ to the DATAJ COMMON block. 
e Delete the statement: 
EQUIVALENCE (T2YJ, IND1), (I2ZJ,IND2). 
The affected subroutines are: CMNGF, CMSET, CMSS, CMSW, 
CMWS, CMWW, EFLD, ETMNS, HINTG, HSFLD, NEFLD, NHFLD, PCINT, 
QDSRC, SFLDS, UNERE. 


1. File 1/O Modifications 
The SOMNTX data file is attached to unit 21. Arrays were written 
to the SOMNTX data file with the statement 
WRITE (21, 100) A 
where A is an array. The array is written correctly to the file, element by element, 
but when NEC-3 attempts to read in the array with the statement 
READ (21, 100) A 
a file read error results. The solution to the problem is to WRITE and READ all 
arrays with an implied DO LOOP. For example, 
WRITE (21, 100) (A(I), I = 1, 10) 
READ (21, 100) (A(1), I = 1, 10). 


The affected subroutine is GNDINO. 
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m. Compiler-specific Subroutines 
Two compiler-specific routines were added to NEC-3. The timer 
subroutine SECOND (ITIME), which returns the INTEGER “4 variable ITIME with 
the number of hundredths of seconds since midnight. The version subroutine 
VERSION prints a line of output which includes the compiler’s name, date, and time 


of the output. Timer and version subroutines are listed in Appendix D. 


2. SOMNTX Modifications 
The following modifications were made to the SOMNTX program, the 


details are the same as for the NEC-3 program: 


e SAVE statements added 

« Generic functions 

¢ FORTRAN 77 standard array passing 
¢ Output statements 

e External functions 

e INTEGER*4 


e File I/O 


It should also be noted that the SOMNTX program was not promoted to 
full double precision for this research, as the accuracy of the output was sufficient 
without doing so. However, the output was extended to double precision prior to 


writing the data file so that the double precision NEC-3 could read the data file. 








IV. RESULTS 


The revised NEC-3 code was evaluated with three types of test problems. A 
set of simple problems, a mixed set of problems from the NEC user’s guide, and a 
set of ground related problems. The NEC-3 input files for all problems tested are 
listed in Appendix E. Since the primary objective of this research was to evaluate 
the feasibility of moving a large mainframe-dependent scientific FORTRAN program 
to a 32-bit personal computer, four versions of NEC-3 (each compiler under a 
different compiler or math co-processor option) were tested for speed and accuracy 
against the IBM 3033AP mainframe computer. 

In addition to the Definicon version, an attempt was made to compile the 
NEC-3 code with all three MS-DOS compilers using both math co-processors. Both 
of the SVS FORTRAN 386 compiler versions, however, repeatedly crashed while 
running the simple 49-segment dipole problem and were not further evaluated. In 
addition, the Lahey compiler does not support the DOUBLE PRECISION 


COMPLEX data type with the Weitek w3167 math co-processor. 


A. DIPOLE EXAMPLES 
Four dipole problems consisting of a 49-segment, a 99-segment, a 199-segment 


and a 299-segment dipole were run on all versions to evaluate the speed and 
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accuracy of each. Table 4.1 compares the full double precision version of NEC-3 
with the single precision version. Tables 4.2-4.5 compare the fill time, factor time, 
and input impedance of each version against the mainframe. Table 4.6 compares the 


execution time of a 299-segment dipole with a 301-segment dipole. 


TABLE 4.1 


SINGLE VERSUS DOUBLE PRECISION 
(LAHEY COMPILER) 


| EXAMPLE | FILL(SEC) | FACTOR INPUT 
(SEC) IMPEDANCE | 
5.66 93 
6.37 ‘99 
21.70 747 
: 24.27 841 
| DIP 199 
| SINGLE 83.93 60.20 
| DOUBLE 94.14 67.45 
| DIP 299 
| SINGLE 186.80 203.50 
| DOUBLE 209.38 228.49 


) TOTAL TIME 


| 
| 
| 
| SINGLE 570.19 | 
| DOUBLE 639.5 | 
| 
| 


| DOUBLE PRECISION TIME PENALTY > = 12% 



















SINGLE 
DOUBLE 


77.927 + j 43.850 
77.900 + j 44.481 










79.041 + j 60.670 
77.997 + ; 44.606 














77.906 + j 39.708 
78.055 + j 44.687 






we 





85.347 + j 12.047 
78.077 + j 44.721 
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TABLE 4.2 


49-SEGMENT DIPOLE 









| SYSTEM FILL TIME | FACTOR INPUT | 
| TIME IMPEDANCE | 


IBM 3033 AP 790 + j4448 | 


i 








AST 386/33 
| LAHEY w/386 77.90 + | 44.48 


1 NDP w/387 7.86 0.77 77.90 + j 44.48 
w/3167 5.93 0.44 77.90 + j 44.48 


DIS-780 77.90 + j 4448 _ | 


1 ——— em ny a ee enn ee 


TABLE 4.3 


99 SEGMENT-DIPOLE 






SYSTEM FILL TIME FACTOR INPUT | 
TIME IMPEDANCE 


IBM 3033 AP 12.25 78.00 + j 44.61 





AST 386/33 | 
LAHEY w/386 2427 | 841s] 78.00 + 5 44.61 | 
NDP w/387 29,93 6.27 78.00 + 544.61 | 

w/3167 22.46 3.9 78.00 + 44.61 | 


DIS-780 117.33 38.22 78.00 + j 44.61 
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TABLE 4.4 
199-SEGMENT DIPOLE 





TABLE 4.5 


299-SEGMENT DIPOLE 


FILL TIME FACTOR INPUT 
TIME IMPEDANCE 


67.36 78.08 + j 44.72 






| IBM 3033 AP 

| 

| 
| LAHEY w/386 209.38 228.49 78.08 + j 44.72 | 
| 


) NDP w/387 255.35 175.93 78.08 + j 44.72 
w/3167 190.04 109.52 78.08 + j 44.72 


| . . : 
| DIS-780 1017.99 1067.64 78.08 + 544.72 | 
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TABLE 4.6 


EXECUTION TIMES OF 299-SEGMENT VS. 301-SEGMENT DIPOLES 
(Lahey Compiler) 


— FILL L TIME | FACT OR TIME | INPUT IMPEDANCE. 
(Q) 


= 209.38 228.49 78.08 + j 44.72 
Leth ee 256.88 418.09 78.08 + j 44.72 


As demonstrated in Table 4.1, the full double precision version of NEC-3 is 









| 





Pome ap tt pee ae re ee eS ERE f 









L. 


only 12% slower than the single precision version. However, a 12% decrease in 
speed is small when compared to the significant improvement in accuracy. Tables 
4.2 through 4.5 show that all double precision versions of the code are equally 
accurate. The two versions using the Intel 80387 math co-processor, Lahey and 
MicroWay, differ in execution speed by less than 1%. The Microway Weitek version 
is 30% faster than the MicroWay Intel 80387 version, with no loss in accuracy. The 
Definicon system is 4.7 times slower than the Intel 80387 versions. 

Table 4.6 illustrated how execution times increase dramatically when the 
number of segments in the problem exceeds the in-core memory limit of 300 
segments. This increase in execution time is directly a function of hard disk speed 


and data transfer rates. While microcomputer hard disk systems have improved 
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dramatically over the past few years, they are still significantly slower than the large 
hard disk systems used by mainframe computers. Because of this, large problems 
requiring an out-of-core solution are best performed on a mainframe computer. 
Although very slow in comparison to the AST 386/33, the Definicon DSI-780 
co-processor board proves that a full double precision 32-bit version of NEC-3 with 
300 segments in-core can be executed in a mere 4 megabytes of RAM. Since the 
Definicon board resides in a host system, the operating system overhead is not 
included in this total. For comparison, assuming the operating system requires a full 
1 megabyte of RAM, the double precision version of NEC-3 with 300 segment in- 
core limit uses only 5 megabytes of RAM on the AST 386/33. If the in-core limit 
was raised to 500 segments, only 7.5 megabytes of RAM would be required, 750 
segments: 12.4 megabytes of RAM, and all 1000 segments would require only 20 
megabytes of RAM. While 20 megabytes of RAM costs a great deal of money, this 


cost is less than half the original cost of the AST computer system. 


B. BASIC ABOVE GROUND PROBLEMS 
The following examples were all taken from the NEC user’s guide, Reference 


13, with the exception of the first example. 
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1. Rhombic Antenna 
This example is a horizontally polarized antenna modeled over a finite 


ground. This significant output data is tabulated in Table 4.7. 


TABLE 4.7 


RHOMBIC ANTENNA HORIZONTALLY POLARIZED 


SEG MENT INPUT IM IMPEDANCE Q 


NUMBER NUMBER Spuumne 


352. Rio faraeeed 06 
352.05 172.06 


HORIZONTAL GAIN NORMALIZATION FACTOR = 17.95 db 





[eres ee eye ee 


t 
| 
| 
| 
| 
\ 
{ 
| 
| 
i 
| 


All versions tested produced exactly the same results. The execution times for 
each problem varied the same between each version as the dipole problems did, thus 


the fill and factor times are not included in the table summaries. 


2. NEC User’s Guide Problems 
The significant results from the first nine problems in the NEC user’s 


guide, Reference 13, are presented here in Tables 4.8 through 4.16. With the 
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exception of problems 3, and 5 the output of all versions of NEC-3 exactly match the 


official validated results listed in the NEC user’s guide. 


TABLE 4.8 


CENTER FED LINEAR ANTENNA 
(Example 1 in Ref. 13) 









| LOAD (Q) INPUT IMPEDANCE (Q) 
REAL IMAGINARY 
CURRENT (A) CHARGE DENSITY 
(C/m) 


PHASE MAG PHASE | 


2.702E-03 -30.68 3.665E-11 60.06 =f 
9.829E-03 -24.34 5.691E-19 24.74 | 


NEAR ELECTRIC FIELDS | 


[MAG (V/m) | PHASE | MAG (V/m) | PHASE _| 





| 1 
| 
=) 




















1.023E-05 1.304E+01 
5.544E +01 1.254E+01 
2.127E +02 4.214E-04 
4.384E +02 3.091E-04 


3.803E +02 








TABLE 4.9 


CENTER FED LINEAR ANTENNA, VARIABLE FREQUENCY 
(Example 2 in Ref. 13) 


NORMALIZED 
IMPEDANCE (Q) 
2.658E +01 - j 6.321E+02 5.31S5E-01 - j 1.264E+01 


4.714E+01 - j 2.724E+02 9.429E-01 - j 5.447E+00 


8.055E+01 + j 4.571E+01 1.611E+00 + j 9.143E-01 


FREQUENCY CONDUCTIVITY INPUT IMPEDANCE 


(MHz) (MHOS/m) 


3.72E+07 1.124E+02 + j 6.543E+01 











TABLE 4.10 


VERTICAL HALF-WAVE DIPOLE OVER GROUND 
(Example 3 in Ref. 13) 


| GROUND | INPUT IMPEDANCE (Q) VERTICAL GAIN 
TYPE NORMALIZATION 
FACTOR (db) 


PERFECT 1.064E +02 + j 9.905E+0 8.52 
FINITE 1.111E+02 + j 1.101E+01 1.54 


NORMALIZED RECEIVING PATTERN 
NORMALIZATION FACTOR = 0.3277E-01 





NOTE = The radiated field near ground for this example differed 
from that listed in Reference 13 by 5% to 10%. 
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TABLE 4.11 


T-ANTENNA ON A BOX OVER PERFECT GROUND 
(Example 4 in Ref. 13) 





TABLE 4.12 


12 ELEMENT LOG PERIODIC ANTENNA IN FREE SPACE 
(Example 5 in Ref. 13) 


INPUT IMPEDANCE (Q) NORMALIZED GAIN FACTOR | 
(db) 


NOTE: Imaginary portion of input impedance is incorrect, 


the correct answer is - j 4.484E-01. 


Radiation patterns differ from that of Reference 13 by 5% to 10%. 
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TABLE 4.13 


CYLINDER WITH ATTACHED WIRE 
(Example 6 in Ref. 13) 


TAG NUMBER - | SEGMENT NUMBER | INPUT IMPEDANCE (a) | 


ae von ase 


ISOLATION DATA 
| 
COUPLING BETWEEN : 
| 
| 















5.493E+01 + j 2527E+01 







SEGMENT SEGMENT 


NUMBER NUMBER NUMBER NU} BER 


| 
a sip tictomil 

| 
| MAXIMUM ieeeecad MAXIMUM COUPLING 
COUPLING 


LOAD IMPEDANCE (22) INPUT 


















(db) 








IMPEDANCE (2) 


1 = -13.709 55.83 - j 20.64 1823 -j 11.64 | 
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TABLE 4.14 


SCATTERING BY A WIRE 
(Example 7 in Ref. 13) 





RADIATION PATTERNS 
FREE SPACE, EXCITATION 0° 
POWER GAIN 


VERT (db) TOTAL (db) _| 


-12.48 -999,99 -12.48 | 
-18.33 -999,99 -18,33 : 


FREE SPACE, EXCITATION 45° 
POWER GAIN 


VERT (db) HORZ (db) 


0 -18.38 - 
45 -10.40 3 


PERFECT GROUND, EXCITATION 45° 


| ANGLE THETA POWER GAIN | 
VERT (db) HORZ (db) TOTAL (db) | 











ANGLE THETA 





> 
a 





i 


ANGLE THETA 






TOTAL (db) | 


-18.38 
-10.40 






S$ 
38 








FINITE GROUND, EXCITATION 45° | 


| ANGLE THETA POWER GAIN : 
VERT (db) HORZ (db) TOTAL (db) | 


NS) 
rs 


















i 

i} 

} oc 
S 
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TABLE 4.15 


STICK MODEL OF AN AIRCRAFT IN FREE SPACE 
(Example 8 in Ref. 13) 


| EXCITATION RADIATION PATTERNS 
: THETA ANGLE POWER GAINS ? 
| ETA 


TH 


VERT (db) | HORZ (db) | TOTAL | 


t 
ei’ 
( 





TABLE 4.16 


BISTATIC SCATTERING BY A SPHERE 
(Example 9 in Ref. 13) 


POWER GAINS 
VERT (db) HORZ (db) | TOTAL an | 


NEAR ELECTRIC FIELDS 
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While the errors in running problems 3 and 5 from the NEC user’s guide 
are significant, they are not catastrophic. In fact, the calculated answers differ from 
the correct answers by only 5% to 10%. What is also unusual about the problem is 
that all versions of the double precision NEC-3 generated exactly the same answers, 
including the incorrect ones. Furthermore, when the mainframe NEC-3 code was 
recompiled with the new IBM VS2 FORTRAN compiler, the results were exactly the 
same as the personal computer versions, including the inaccuracies of Problems 3, 
and §. This demonstrates that the preblem is not hardware or compiler related, but 
rather related to how NEC-3 was written. It is possible the new IBM (VS2) and 
personal computer FORTRAN 77 compilers somehow implement data structures 
such as COMMON blocks and EQUIVALENCE statements in a slightly different 
manner than the old IBM (VS1) compiler. It is also possible that the problem is a 
result of a non-rigorously written section of code. Since much of the NEC-3 program 
was written in the style of FORTRAN 66, the less tolerant FORTRAN 77 compilers 


might have incorrectly compiled a small section of the NEC-3 code. 


C. GROUND PROBLEMS 
Tables 4.17 through 4.21 tabulate the significant results of several ground 
related problems. As with the previous example, all versions of NEC-3 produced 


exactly the same correct results. 
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TABLE 4.17 





MONOPOLE OVER A RADIAL WIRE SCREEN OVER FINITE GROUND 







| SCREEN HEIGHT | INPUT IMPEDANCE | AVERAGE POWER | 
| (m) (Q) GAIN | 
| 0.2 83.05 - j 3346.5 0.5468 : 
| 0.6 224.1 - j 86.80 0.3876 | 





53 





TABLE 4.18 
DIPOLE ANTENNA MOVING CLOSER TO GROUND 


HEIGHT ABOVE INPUT IMPEDANCE (Q) 


GROUND (m) 


0.1 66.34 + j 62.50 
0.03 73.89 + j 75.49 
0.01 97.87 + j 127.95 
0.003 119.87 + j 219.62 
0.001 132.87 + j 322.13 
0.0003 144.21 + j 441.81 
0.0001 154.59 + j 554.56 
0.00003 167.34 + j 681.82 
0.00001 180.94 + j 802.05 
3.0E-6 198.95 + j 939.39 
1.0E-6 219.38 + j 1071.00 
3.0E-7 249.11 + j 1223.50 *** 


1.0E-7 280.04 + j 1361.90 *** 


*** NOTE: Integration was step size limited, although impedance 


calculations are correct. 
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TABLE 4.19 


DIPOLE IN DIELECTRIC AND LOSSY MEDIA 


CONDUCTIVITY = 0.0 
INPUT IMPEDANCE = 11.50 - j 32.66 


RADIATION PATTERNS 


| ANGLE THETA POWER GAINS 
VERT (db) HORZ (db) TOTAL (db) 


| 
| 
| 


CONDUCTIVITY = 0.89E-02 
INPUT IMPEDANCE = 70.34 + j 11.02 


RADIATION PATTERNS 


| ANGLE THETA POWER GAINS 
VERT (db) HORZ (db) TOTAL (db) 


| 
| 
| 
| 
| 
| 
| 
| 
1 
| 
| 
| 
| 
| 
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TABLE 4.20 


MONOPOLE ON A GROUND STAKE 


INPUT IMPEDANCE = 96. 14 +j 38 Al (Q) 
AVERAGE POWER GAIN = 0.3215 db 


NEAR ELECTRIC FIELDS 


PHASE | MAG(V/m)| PHASE 


0.172E-05 0.106E-04 
0.121E-05 : 0.112E-04 
0.116E-05 : 0.121E-04 
0.103E-05 ; 0.223E-04 
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TABLE 4.21 


MONOPOLE ON A 6-WIRE RADIAL GROUND SCREEN 


INPUT IMPEDANCE = 53. 96 + j j 29. 45 (Q) 


AVERAGE POWER GAIN = 0.5821 (db) 


0.293E-05 0.179E-04 
0.282E-05 . 0.163E-04 
0.178E-05 : 0.258E-04 
0.174E-05 : 0.276E-04 


0.175E-05 ; 0.380E-04 
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Vv. CONCLUSIONS 


The results of testing the double precision version of NEC-3 with a wide 
variety of problems, clearly show that NEC-3 can be successfully implemented on 
a 32-bit personal computer. 

The AST Premium 386/33 utilizing the Intel 80387 math co-processor is 2.3 
times slower than the IBM 3033 mainframe. When using the Weitek w3167, the AST 
Premium 386/33 is only 1.62 times slower than the mainframe. Since the fill and 
factor times of the IBM 3033AP mainframe are true processor times, without the 
obvious time-sharing overhead added to them, the speed of the AST Premium 
386/33 is very close to the practical speed of a moderately loaded mainframe 
computer. 

The full double precision NEC-3 code exactly matched the output of the IBM 
3033AP mainframe for accuracy. The additional time required to perform all 
calculations in double precision arithmetic was shown to be only 12% more than the 
full single precision version. 

The Definicon DSI-780 demonstrated that the full double precision NEC-3 can 
be successfully executed in 4 MB of RAM. When the operating system is taken 


into account, the full double precision NEC-3, with 300 unknowns in-core, can 
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be operated on any 32-bit personal computer with at least 5 MB of RAM. If the in- 
core limit were raised from the current 300 unknown to 1000 unknowns, the double 
precision NEC-3 code could be operated on a system with at least 20 MB of RAM. 
Due to the inherent slow speed of microcomputer hard disk systems relative to the 
mainframe, it is recommended that the NEC-3 code be converted to the largest 
possible number of unknowns in-core for use on personal computers. 

The inaccuracies discovered while running the sample problems from the NEC 
User’s Guide indicate that the “bug" is not hardware or software dependent, but a 
matter of how the new FORTRAN compilers are implemented. It is possible that 
the solution to the problem may lie in how the new FORTRAN compilers 
implement shared memory, since NEC-3 uses shared arrays via COMMON blocks 
and EQUIVALENCE statements throughout the program. It is recommended that 
the code be rewritten to eliminate the use of EQUIVALENCE statements for 


reusing arrays. 
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APPENDIX A 


LAHEY F77L3 FORTRAN COMPILER 


Switches and Options 


IMPLICIT NONE in effect 

Standard implicit typing in effect 

Load NDP (80287 or 80387) nonINTEGER operands 
from memory 

NDP nonINTEGER operands may be kept in the 
NDP between statements 

Bounds check array subscripts 

No bounds checking 

Output warning messages indicating statements that 
do not conform to the FORTRAN 77 Standard 
Suppress the output of warning messages indicating 
nonconformances to the FORTRAN 77 Standard 
Direct files are processed without headers 

Direct files require F77L-EM/32-generated headers 
Free-format source file 

Standard-format source file 

Hardcopy (source listing) output 

No hardcopy (listing) output 

Subprogram interface check 

No subprogram interface check 

Generate Weitek 1167 and 3167 code 

No 1167 and 3167 code generated 

Line-number table 

No line-number table 

Output compiler options 

No options output 

Protect constant arguments 

No protection 

Ensure code never uses more that 8 NDP registers at 
a time 

Infinite NDP stack for expressions 

Create the protected-mode side of an RPC interface 
module 











Don’t create the protected-mode side of an RPC 
interface module 

Create the real-mode side of an RPC interface 
module 

Don’t create the real-mode side of an RPC interface 
module 

Remember (SAVE) subprogram local variables and 
arrays. 

No remembering local items 

Generate SOLD3 information file 

No SOLD3 information file generated 

Type; default length of INTEGER type is 2, 
LOGICAL type is 1 

Default length of INTEGER is 4, LOGICAL is 4 
Warning messages output 

Not all warning messages output 

Cross-reference listing 

No cross-reference listing 

Perform production code optimizations (limits 
SOLD3 and P77L Profiler interfaces) 

No optimizations that would limit SOLD3 and P77L 
Profiler interfaces 
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APPENDIX B 


MICROWAY NDP FORTRAN-386 COMPILER 


-c 
-list 
-LIST 


-Idir 
-o filename 


-R 


-S 


“Vv 


Switches and Options 


Compile to the obj” level only 

Create a ".ist" file with errors annotated 

Same as above, but listing does not show 
pathnames for source 

Alternate directory(ies) for INCLUDE files 
Place executable file output into the file named: 
"filename" 

Put all data in the code segment, i.e., generate 
"ROMable" code 

Do not produce object files or executable files 
Have the compiler driver print out the program 
name and command-line arguments as it runs each 
subprocess 

Suppress warning messages 

Search the library named c:\NDP\libname.lib 
Place the symbols command into the linker file for 
Phar Lap linker 

Turn on trace/debug information 

Generate a frame pointer for stack traces 
Allow arbitrary file names to be specified to the 
compiler 

Display the names of files as they are opened 
Do not stop in the event of a code generator abort 
Do not use 386 inline library routines 

80287 inline code with library cails for 
transcendentals 

80287 inline code with inline transcendentals 
(except sin and cos) 

80387 inline code with all transcendentals 
(including sin and cos) 

Advanced 80387 stack utilization 

Weitek W1167 code 

Produce W1167 macro instructions 
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-n6o 
-CG1 


CG2 
-CG3 
-CG4 
-CG5 


-CG6 

-CG7 

-CG8 

-CG9 

-CG10 

-O 

-OL 

-OM or -OLM 


-ONW 
-OFF 


-OFFP 
-OFFR 


-OFFA 
-OFFH 
-OFFS 


-OFFN 
-i2 
-onetrip 
“u 

-ep 

-U 


-fl 





Promote no-Float switch for Weitek 

Turn on runtime checking of subranges and array 
bounds 

Allocate all variable to memory 

Allocate code temporaries to memory 

Prepend all variables with an underscore 

Output an assembly file using ".asm" extension; the 
default is ".s" 

Do not put an underscore in front of the names of 
global variables and procedures 

Avoid jumps with inline code 

Turn of jump shortening 

Lengthen bytes in jump 

Suppress version number in data segment 

Perform default optimizations 

Add speed optimizations related to moving code 
out of loops and speeding up loops in general 
Add additional memory optimizations to -O and - 
OL 

Emit a warning when dead code is eliminated 
Activates -OFFP, -OFFR, -OFFA, -OFFH, -OFFS 
and -OFFN 

Disable peephole optimizer 

Do not allocate programmer-defined local variable 
to a register unless they are declared register 

Do not move frequently used procedure and data 
addresses into registers 

Turn off the code-hoisting optimization 

Turn off the optimization that deletes all code that 
stores into or modifies variables that are never 
read from 

Do not move invariant floating-point expressions 
out of loops 

Make the type INTEGER be INTEGER*2 
Execute at least one iteration of every DO loop 
Make "undefined’ the default data type for 
undeclared variables 

Allow inline comments with "!" as in VAX 
FORTRAN 

Do not convert upper-case names in FORTRAN to 
lower-case 

This makes characters unsigned as they are in 
some implementations of FORTRAN 
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Turn of compiler-time checking of FORMAT 
statements 

Pad Hollerith constants on the right with blanks 
Compile lines starting with "x", "X", "d", "D". 
Do not allow dollar signs in names 











+q -q 
+p -p 
+X 

+d -d 
+c72 


+save 


+saveall 
+charequ 


+int2 


+log2 
+de -de 


+f -f 
+287 
+387 
+w1167 
-lfname 
-e fname 


-ifname 


APPENDIX C 


SVS FORTRAN COMPILERS 


Switches and Options 


Show more (-q) or less (+q) information on the 
progress of the compile to the user 

Prompt (+p) or don’t prompt (-p) to the standard 
input in the case of a compile time error 
Generate a cross reference in the listing file 
Controls whether the compiler does (+d) or does 
not (-d) generate tables of information for the SVS 
Debugger 

Truncate input lines to 72 columns 

Compile program as if each subprogram contained 
a SAVE statement 

Same as setting the SSAVEALL option 

Allow certain non ANSI Standard associations 
between character and numeric data 

Change default attribute of INTEGERs and 
LOGICALs 

Change default length attribute of LOGICALs 
Conditionally compile debugging source code 
found in lines with the letter D in column 1. 
Generate code for floating point hardware (+f) or 
floating point software (-f) 

Generate Intel 287 floating point code, on Intel 
386 only 

Generate Inte] 387 floating point code, on Intel 
386 only 

General Weitek 1167 coprocessor floating point 
code, on Intel 386 only 

Create a listing file of the source program in the 
file named fname 

Place a summary of the compile time errors on file 
named fname 

Name the ".i" file fname 
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+b -b 


+e -¢€ 


+cc -CC 
+a -a 


+c -C 


Control interpretation of the "CUNFORMATTED’ 
file FORM in OPEN statements 

The +e and -e option control whether the 
arguments to subroutines and functions are popped 
in the exit code of the routine before return (+e) 
or by the caller after return from the routine (-e) 
Control of parameter passing order and poping of 
parameters in subroutine calls 

Flag non-ANSI Standard features which are 
utilized by the program 

Setting +c instructs the compiler to create an 
information file of the source file positions of 
compile time errors 











APPENDIX D 


TIMER AND VERSION SUBROUTINES 


1. LAHEY F77L3-EM/32 
INCLUDE C:\NDP\DNEC3.F 


SUBROUTINE VERSION 

CHARACTER*11 CTIME 

CHARACTER*8 CDATE 

CALL DATE(CDATE) 

CALL TIME(CTIME) 

WRITE (6,100) CDATE,CTIME 

100 FORMAT (//,36X,’VERSION - LAHEY DOUBLE W/SAVES  DATE:',A,3X 

&,'TIME:',A,//) 

END 


SUBROUTINE SECOND( TIME) 
INTEGER*4 TIME 

CALL TIMER( TIME) 
RETURN 

END 








2. MICROWAY NDP FORTRAN386 


INCLUDE "DNEC3.F" 


SUBROUTINE VERSION 

CHARACTER*11 CTIME 

CHARACTER*8 CDATE 

CALL DATE(CDATE) 

CALL TIME(CTIME) 

WRITE (6,100) CDATE,CTIME 

100 FORMAT (//,36X,’VERSION - MICROWAY DOUBLE W/SAVES DATE:’,A,3X 

&,‘'TIME:',A,//) 

END 


SUBROUTINE SECOND( ITIME) 

INTEGER*4 ITIME,YEAR,MONTH, DAY,HOUR , MINUTE , SECOND ,MICROSEC 
CALL TIMEDATE( YEAR , MONTH, DAY , HOUR ,MINUTE, SECOND , MICROSEC) 
ITIME = HOUR*360000 + MINUTE*6000 + SECOND*100 + MICROSEC*1E-4 
RETURN 

END 








3. SVS FORTRAN386 


S$INCLUDE C:\NDP\DNEC3.F 


100 


SUBROUTINE VERSION 

WRITE (6,100) 

FORMAT (//,36X,‘VERSION - SVS (MSDOS) W/SAVES ‘,//) 
END 


SUBROUTINE SECOND(ITIME) 

INTEGER*4 ITIME, IHR, IMIN, ISEC,1100 

CALL GETTIM(IHR, IMIN, ISEC,1100) 

ITIME = IHR*360000 + IMIN*6000 + ISEC*100 + 1100 
RETURN 

END 
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4. SVS DEFINICON FORTRAN 


S$SYSTEM 
SINCLUDE C:\DNEC3.F 


100 


SUBROUTINE SECOND(ITIME) 
INTEGER*4 ITIME,HOUR,MIN,SEC,HUND, TIME 


TIME= STIME() 

HOUR = ISHFT(TIM:., -24) 

MIN = ISHFT(IAND(16711680, TIME) ,-16) 

SEC = ISHFT(IAND(65280, TIME), -8) 

HUND = IAND(255,TIME) 

ITIME = HOUR*360000 + MIN*6000 + SEC*100 + HUND 
RETURN 

END 


SUBROUTINE VERSION 
WRITE (6,100) 
FORMAT (//,36X, ‘VERSION - SVS-DEFINICON DOUBLE W/SAVEALL 


END 
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APPENDIX E 


NEC TEST PROBLEMS 


DIPOLE WITH 49 SEGMENTS 
1,49,0,0,0,0,0,.5, .00001 


0,1,25 
-1 


DIPOLE WITH 99 SEGMENTS 
1,99,0,0,0,0,0,.5, .00001 


0,1,50 
-1 


DIPOLE WITH 199 SEGMENTS 
1,199,0,0,0,0,0,.5, .00001 


0,1,100 
-1 


DIPOLE WITH 299 SEGMENTS 
1,299,0,0,0,0,0,.5, .00001 


0,1,150 
-1 
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CM RHOMBIC ANTENNA HORIZONTALLY POLARIZED 

CM LEG LENGTH=398.0 FT. 

CM CENTER WIDTH=314.0 FT. 

CM APEX ANGLE=44.0 DEGREES. 

CM HEIGHT ABOVE GROUND=160.0 FT. 

CM GROUND PARAMETERS-EPSILON=80. SIGMA=4. MHOS/M. (SEA WATER) 
CE CONDUCTOR-AWG NO. 10 WIRE DIA.=0.00425 FT. 
GW1,40,0.0,0.0,160.0,366.082,157.0,160.0,0.00425 

GW2 ,40, 366.082 ,157.0,160.0,732.164,0.0,160.0,0.00425 


GX2,010 
GS0,0,0.304801 

GE1 

FRO,0,0,0,10.0 
GNO,0,0,0,80.0,4.0 
LDc ,2,40,40,300.0 
LDO,4,40,40,300.0 
EX0,1,1,0,0.5 


EX0,3,1,10,-0.5 
RPO,37,1,1401,90.0,0.0,-5.0,0.0 
EN 


TZ 





CEEXAMPLE 1. 


Gwo, 7, 0., 0., -.25, 0., 


GE 

EXO, 0, 4, QO, 1. 
XQ 

LpO, 0, 4, 4, 10., 
PQ 


NEO, 1, 1, 15, .001, 0., O. He 

NX 

CMEXAMPLE 2. CENTER FED LINEAR ANTENNA. 

CM CURRENT SLOPE DISCONTINUITY SOURCE. 
CM 1. THIN PERFECTLY CONDUCTING WIRE 
CE 2. THIN ALUMINUM WIRE 
GWO,8,0.,0.,-.25,0.,0.,.25, .00001 

GE 

FRO,3,0,0,200. ,50 

EX5,0,5,1,1.,0.,50 

XQ 

LD5,0,0,0,3.720E+07 

FRO,1,0,0,300 

EX5,0,5,0,1 

XQ 

NX 

CMEXAMPLE 3. VERTICAL HALF WAVELENGTH OVER GROUND 
CM EXTENDED THIN WIRE KERNEL USED 

CM 1. PERFECT GROUND. 

CM 2. IMPERFECT GROUND INCLUDING GROUND WAVE AND RECEIVING 
CE PATTERN CALCULATIONS 
GWO,9,0.,0.,2.,0.,0.,7.,.3 

GE1 

EK 

FRO,1,0,0,30. 

EX0,0.5,0,1. 

GN1 


RPO,10,2.1301,0.,0.,10. ,90. 


GNO,0,0,0,6.,1.000E-03 


RPO,10,2,1301,0.,0.,10.,90. 
RP1,10,1,0,1.,0.,2.,0.,10.000E+04 
.10. 


EX1,10,1,0,0.,0.,0. 
PT2,0,5,5 

XQ 

EN 


3.E-09, 5.3E-11 





CENTER FED LINEAR ANTENNA. 
.001 
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CEEXAMPLE 4. T ANTENNA ON A BOX OVER PERFECT GROUND 
SPO,0,.1,.05,.05,0.,0.,.01 
SPO,0,.05,.1,.05,0.,90.,.01 

Gx0,110 
$p0,0,0.,0., 
GW1,4,0.,0., 
GW2,2,0.,0.,. 
GW3,2,0.,0., 
GE1l 

GN1 
EX0,1,1,0,1. 
RPO,10,4,1001,0.,0.,10.,30. 


45. o-. "3, .001 
1 


CM EXAMPLE 5. 
CM 12 ELEMENT LOG PERIODIC ANTENNA IN FREE SPACE. 


CM 78 SEGMENTS. SIGMA=D/L RECEIVING AND TRANS. PATTERNS. 


CM DIPOLE LENGTH TO DIAMETER RATIO=150. 


CE TAU=0.93. SIGMA=0.70. BOOM IMPEDANCE=50. OHMS. 


GW 1,5,0.,-1.,0.,0.,1.,0.,0.00667 
GW 2,5,-.7527,-1.0753,0.,-.7527,1.0753,0., .00717 


GW 3,5,-1.562,-1.1562,0.,-1.562,1.1562,0., .00771 
GW 4,5,-2.4323,-1.2432,0.,-2.4323,1.2432,0., .00829 
GW 5,5,-3.368,-1.3368,0. ,-3.368,1.33368,0., .00891 
GW 6,7,-4.3742,-1.4374,0. ,-4.3742,1.4374,0., .00958 
GW 7,7,-5.4562,-1.5456,0. ,-5.4562,1.5456,0., .0103 
GW 8,7,-6.6195,-1.6619,0. ,-6.6195,1.6619,0., .01108 
GW 9, 7,-7.8705,-1.787,0.,-7.8705,1.787,0., .01191 


GW10,7,-9.2156,-1.9215,0. ,-9.2156,1.9215,0., .01281 
GW11,9,-10.6619,-2.0662,0. ,-10.6619,2.0662,0.,.01377 
GW12,9,-12.2171,-2.2217,0. ,-12.2171,2.2217,0., .01481 


GE 

FR 0,0,0,0,46.29 

TL 1,3,2,3,-50. 

PL. 2595999 6200. 

TL 3,3,4,3,-50. 

TL 4,3,5,3,-50. 

TL 5,3,6,4,-50. 

TL 6,4,7,4,-50. 

TL 7,4,8,4,-50. 

TL 8,4,9,4,-50. 

TL 9,4,10,4,-50 
TL10,4,11,5,-50. 
TL11,5,12,5,-50.,0.,0.,0.,.02 
EX 0,1..3,20; 17. 

RP 0,37,1,1110,90.,0.,-5.,0. 
EN 
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CM EXAMPLE 6 
CE CYLINDER WITH ATTACHED WIRES. 


$P0,0,10.,0.,7.3333,0. 
$P0,0,10.,0.,0.,0.,0., 


,0.,38.4 
38.4 


SPO,0,10.,0.,-7.3333,0.,0.,38.4 


GM0,1,0.,0.,30. 
S$P0,0,6.89,0.,11.,90., 


0. ,44.88 


SPO,0,6.89,0.,-11.,-90.,0. 44.88 


GRO,6 
SP0,0,0.,0.,11.,90.,0. 


44.89 


S$P0,0,0.,0.,-11,-90.,0. ,44.89 


GW1,4,0.,0.,11.,0.,0., 


23.,.1 


GW2,5,10.,0.,0.,27.6,0.,0.,.2 


GS0,0,.01 
GE? 
FRO,1,0 
CP1,1.2 
1 


= | 
1000,0.,0.,5. 
1 


CM SAMPLE PROBLEM NO. 
CM 1. STRAIGHT WIRE - 
CM 2. STRAIGHT WIRE - 
CM 3. STRAIGHT WIRE - 
CE (SIGMA=.0001 
GWO,15,-55.,0.,10.,55. 
GE1l 

FRO,1,0,0,3. 
EX1,2,1,0,0.,0.,0.,45. 
RPO,2,1,1000,0.,0. ,45. 
GN1 
EX1,1,1,0,45.,0.,0. 


7 
FREE SPACE 





PERFECT GROUND 


FINITELY CONDUCTING GROUND 


EPSILON=6 ) 
;0.,10;., 01 
0. 

10% 


RPO,19,1,1000,90.,0.,-10.,0. 


GNO,0,0,0,6.,.0001 


RPO,19,1,1000,90.,0.,-10.,0. 


EN 
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CM SAMPLE PROBLEM NO.8 
CE STICK MODEL OF AN AIRCRAFT IN FREE SPACE 
GW1,1,0.,0.,0.,6.,0.,0.,1. 

GW2,6,6.,0.,0.,44.,0.,0.,1. 


CM BISTATIC SCATTERING BY A SPHERE 
CM PATCH DATA ARE INPUT FOR A SPHERE OF 1. METER RADIUS 
CM THE SPHERE IS THEN SCALED SO THAT KA=FREQ. IN MHZ. 

CM THE PATCH MODEL MAY BE USED FOR KA LESS THAN ABOUT 3. 
CE FOR THIS RUN, *** 
SPO0,0, .13795, .13795,. 
SPO,0, .51328, .21261,. 
.51328,. 
.21520,. 
.58794,. 
.80314,. 
.19134,. 
.54490,. 
.81549,. 
.96194,. 


SPO,0,.21261, 
SPO,0, .80314, 
SPO,0, .58794, 
SPO0,0,.21520, 
SP0,0, .96194, 
SPO,0, .81549, 
SPO,0, .54490, 
SP0,0, .19134, 
Gx0,111 


GS0O,0,47.71465 


GE 
FRO,1,0,0,2.9 


EX1,1,1,0,90.,0.,0. 


KA=2 .9 . *** 


98079,78. 
83147,56. 
83147, 56. 
55557,33. 
55557,33. 
95957,,33.. 
19509,11. 


19509,11 
19509 11 
19509,11 


RPO ,19,1,1000,90.,0.,-10.,0. 
RPO,1,19,1000,90. ,0.,0.,10. 


NEO,1,1,11,0. 
NEO,1,11,1,0. 
NEO,11,1,1,0. 
NHO,1,1,11,0. 
NHO,1,11,1,0. 
NHO,11,1,1,0. 
EN 


ooo0oo0o0°o 


’ ? 
’ ’ 
9 oar} 
’ » 
’ 


owoouwo 
oowoowm 


75,45. 
25,22. 


25,67 


75,15. 
75,45. 
75,75. 
257,11. 
<2Dy oo 
(255.90: 
.25,78. 


,.11975 
5,.17025 
.3,.17025 

, 16987 

,- 16987 

,- 16987 
25, .15028 
75, .15028? 
25, .15028 
75, .15028 
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CM TEST OF THE NEW NEC WITH SOMMERFELD 


CM GREEN'S FUNCTION FOR RADIAL WIRE SCREEN OVER FINITE GROUND 
CM SCREEN RADIUS=60. M (ONE WAVELENGTH RADIUS) 
CM SCREEN HEIGHT=.02M 6 RADIAL WIRES 


GRO, 6 


CM MONOPOLE ON RADIAL WIRE GROUND SCREEN 


GW1,6,0.,0.,.01,0.,0.,15.2, .003 
GEl 

EX0,1,1,0,1. 
RPO,19,2,1001,0.,0.,5.,90. 

NX 

CM 

cM 

CM NOW THE SCREEN JACKED UP TO .6M HEIGHT 
CE 
GWO,12,0.,0.,.6,30.,0.,.6, .003 
GRO, 6 


CE MONOPOLE ON RADIAL GROUND SCREEN 
GF 

GW1,6,0.,0.,.6,0.,0.,15.6, .03 

GE1l 

Ex0,1,1,0.1. 
RPO,19,2,1001,0.,0.,5.,90 

EN 








CM TEST OF NEC GROUND/ DIPOLE CLOSE TO GND, MOVING EVER CLOSER 
CE 
GWO,5,-.25,0.,.1,.25,0.,.1,1.E-10 


GE 

GN2,0,0,0,4. ,-1.7976 
Ex0,0,3,0,1. 

XQ 

NX 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.03,.25,0., .03,1.E-10 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.01,.25,0.,.01,1.E-10 


NX 
CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.003,.25,0., .003,1.E-10 


NX 

CE TEST OF NEC GROUND 
GWO.5,-.25,0.,.001,.25,0.,.001,1.E-10 
GE 


CE TEST OF NEC GROUND 
GWO,5,-.25,0., .0003,.25,0., .0003,1.E-10 
GE 


NX 

CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.0001,.25,0.,.0001,1.E-10 
GE 
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CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.00003, .25,0., .00003,1.E-10 


GN2,0,0,0,4.,-1.7976 
0,3,0, 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,.00001, .25,0.,.00001,1.E-10 


GN2,0,0,0,4. ,-1.7976 
Ex0,0,3,0,1 

XQ 

NX 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,3.E-6,.25,0.,3.E-6,1.E-10 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,1.E-6,.25,0.,1.E-6,1.E-10 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,3.E-7,.25,0.,3.E-7,1.E-10 


CE TEST OF NEC GROUND 
GWO,5,-.25,0.,1.E-7,.25,0.,1.E-7,1.E-10 
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TEST OF NEC3 GROUND/BURIED WIRES & SUCH STUFF 
1. DIPOLE IN DIELECTRIC AND LOSSY MEDIA 
TEST OF UM CARD 


1,7,0,0,-1.5,0,0,1.5,.01 


,0,0,16 
0,1,1000,0,0,10,0 
,0,0,16,-16 
0,1,1000,0,0,10,0 


2. MONOPOLE ANTENNA ON A GROUND STAKE 


1,8,0,0,-2.,0,0,0,.01 
2,10,0,0,0,0,0,15, .01 


9,2,1001,0,0,5,90 
/1,21,5000,0,0.0,0,10 


3. MONOPOLE ON 6-WIRE RADIAL GND SCRN 
FIRST THE WIRE GROUND SCREEN 


14,12.0,-.05,.8,0,-.05, .01 
1,.8,0,-.05,0,0,0,.01 
6 


NOW THE MONOPOLE 
2,10,0,0,0,0,0,15, .01 
1 


0, 2% 
0,19,2,1001,0,0,5,90 
0,1,1,21,5000,0,0,0,0,10 


11. 


12. 
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